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We describe a model independent search for physics beyond the standard model in lepton final 
states. We examine 117 final states using 1.1 fb _1 of pp collisions data at yfs = 1.96 TeV collected 
with the DO detector. We conclude that all observed discrepancies between data and model can be 
attributed to uncertainties in the standard model background modeling, and hence we do not see 
any evidence for physics beyond the standard model. 
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I. INTRODUCTION 

The standard model (SM) has been remarkably suc- 
cessful in accommodating all the interactions between the 
fundamental particles Despite this success, there are 
strong motivations to expect new phenomena at ener- 
gies at the order of the electroweak scale. For example, 
the Higgs boson [2[ receives quantum corrections to its 
mass through loop diagrams. The scalar nature of the 
Higgs boson leads to a quadratic divergence, with an up- 
per limit of the integral set by the highest scale, i.e., the 
Planck mass (10 19 GeV). To maintain the Higgs mass 
close to the electroweak scale, it is necessary to fine tunc a 
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jmrameter in the theory to within Mi^/A/pianck ~ 10 16 

There are few logical options for overcoming this prob- 
lem. If the Higgs boson does not exist, then there must 
be a new contribution to the physics at the electroweak 
scale. If the Higgs boson docs exist, then the theory 
must be either fine tuned or a generalized Higgs scheme, 
beyond the SM, is present at the electroweak scale. 

Assuming that beyond standard model (BSM) physics 
exists, we do not know how it appears, rendering its 
search difficult. While there are many theories that 
predict observable differences with the SM, these mod- 
els usually depend on additional unspecified parameters 
which broaden the possible range of results. 

Motivated by uncertainty and expectations of physics 
beyond the SM, we examined data from many channels 
in pp collisions at y/s = 1.96 TcV at the Tevatron Col- 
lider at Fcrmilab, collected by the DO experiment, for 
deviations from the SM. After this, we focus on events 
with objects with high transverse momentum (j>t) in a 
quasi-model-independent search of new phenomena ef- 
fects. Our background model is specific for final states 
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containing lcptons, which form the focus of this paper. 
Similar approaches have been applied to data from the 
DO Collaboration the HI Collaboration at the 

HERA ep collider at DESY 0, and the CDF Collab- 
oration at the Tevatron Q . 

Our technique trades the sensitivity of specific searches 
for breadth of coverage: we do not design selections fo- 
cused on a particular model and neglect systematic un- 
certainties. This way, we can incorporate many channels 
without developing a detailed modeling for each individ- 
ual channel. This approach limits sensitivity for physics 
beyond the SM in individual final states, but it helps 
identify global differences relative to the SM expecta- 
tions. If any particular final state or distribution found 
discrepant with the SM remains significantly discrepant 
after systematic uncertainties are considered, then it war- 
rants claim for the presence of physics beyond the SM. 
The benefit of this approach is that we can look in a co- 
ordinated way at many channels, applying expectations 
from the SM and a model of the detector in a relatively 
straightforward manner, to search for discrepancies be- 
tween data and the SM. 

The data for the search consists of events containing 
objects that have large px values. We divide the data 
and the selected Monte Carlo (MC) simulated events into 
seven inclusive subsets based on the number and types of 
lcptons identified in each event. Unlike the search con- 
ducted by the CDF Collaboration d, Q , only events with 
at least one electron or muon are considered. For each of 
the chosen final states, we apply corrections to the MC 
simulation, as determined from the previous DO studies 
based on well-understood regions of phase space, domi- 
nated by particular SM processes, as discussed in Sec.lIVI 
To account for any incorrect normalizations in the ab- 
sence of systematic uncertainties, we fit for contributions 
from each of the subsets to obtain scale factors which re- 
produce the distributions in the selected data with MC 
events and multijet background events determined from 
data, as discussed in Sec. |V] 

The seven non-overlapping inclusive subsets are 
merged to create an input file for the analyses employing 
algorithms called vista and sleuth [8( , as discussed in 
Sec. ED 



II. DO DETECTOR 

The data correspond to 1.07 ± 0.07 fb _1 of integrated 
luminosity from pp collisions at the Tevatron Collider at 
Fermilab, collected with the DO detector at y 7 ! = 1.96 
TeV during 2002-2006. 

The DO detector is described in detail elsewhere fl(ij |. 
The central tracking, calorimctry, and muon systems are 
the components most important to this analysis. The 
central tracking system consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both 
located within a 2 T superconducting solenoidal magnet, 
and provides charged particle tracking for pscudorapidi- 



ties I77I < 3, where r\ = — ln[tan(6>/2)], and 9 is the polar 
angle relative to the center of the detector with respect 
to the proton beam direction. 

The three liquid-argon/uranium calorimeters are 
housed in separate cryostats. Outside of the tracking 
system, a central section covers up to \rj\ = 1.1. Two 
end calorimeters extend coverage to = 4.2. The 
calorimeter is highly segmented with four electromag- 
netic (EM) and four to five hadronic longitudinal layers; 
transverse to the particle direction, typical segmentation 
is Ai] = A<p = 0.1, where (f> is the azimuthal angle. 

Beyond the calorimeter, a muon system consists of 
a layer of tracking detectors and scintillation trigger 
counters in front of 1.8 T iron toroids, followed by two 
similar layers after the toroids, all at pseudorapiditics 



< 2.0 



11] 



A three level trigger system selects events, recording 
data at about 100 Hz. Our sample was collected using 
triggers that select events with at least one electron or 
one muon. 



III. OBJECT ID AND EVENT SELECTION 

In this section, we describe the identification criteria 
used to select energetic objects isolated from other event 
activity, viz., electrons (e ± ), muons (z^), tau leptons 
(t ± ), missing transverse energy ($t), jets, and 6-quark 
jets. The selection criteria for all these objects are identi- 
cal for all final states. In addition, we discuss the criteria 
for non- isolated electrons and muons, i.e., objects that 
are not truly isolated but can mimic the signatures of 
isolated leptons. Because of the difficulty of modeling 
such false leptons, their contributions are estimated di- 
rectly from data. 



A. Vertices 

Only pp interaction vertices reconstructed from at least 
three tracks are allowed in this analysis. Based on the 
Pt of the tracks associated with that vertex, we define 
the primary pp interaction vertex (PV), as the one with 
smallest probability of originating from a minimum-bias 
interaction pjj]. The z coordinate of the PV (zpv) is 
required to be |zpy| < 60 cm (where the z axis is the 
axis along beam direction, with origin at the center of 
the detector). 



B. Electrons and Photons 

Electrons are characterized by an isolated shower in the 
calorimeter and an isolated track in the central tracker. 
Starting with a seed cell, a calorimeter cluster is formed 
using cells within a cone of radius A7Z < 0.4 where 
A1Z = -J {At]) 2 + (A(f)) 2 . Such clusters are required 
to pass the calorimeter isolation criterion (E to t(AlZ < 
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0.4) - £ EM (A7£ < 0.2))/ J B EM (A^ < 0.2) < 0.2, where 
E tot is the total energy of the shower, summing the EM 
and hadronic calorimeter cells, and -Eem is the energy in 
the EM calorimeter only. Every accepted cluster must 
have 90% of E to t within the EM calorimeter, pass a x 2 - 
based selection on the spatial distribution of the shower, 
and be matched with a track extrapolated from the cen- 
tral tracker. An electron likelihood (L e ), based on seven 
tracking and calorimctric parameters, is used to enhance 
signal purity of the candidate electrons. Photons are 
identified as electromagnetic clusters that pass the same 
isolation and shower criteria, but fail to match with a 
track. 

In this analysis, we use only electrons that are found in 
the central calorimeter (CC), with \r)\ < 1.1 and pr > 15 
GeV. Typical electron detection efficiencies are 70% to 
80%. 

To estimate the contribution from non-isolated elec- 
trons (e.g. from multijet background), we use the same 
selection as for signal, but with a reversed L e likelihood 
criterion. 



decays r — > nv T are referred to as Type-1. Decays corre- 
sponding to t~^~ — > 7r^~ titt^ v T are referred to as Type-2 (n 
is an integer > 1), and decays to multiple charged pions 
are referred to as Type-3 decays. Type-3 decays differ 
from Type-1 (ti) and Type-2 (T2) by being matched to 
multiple tracks, and are not used in this analysis. Type- 
1 and Type-2 decays are required to have |?7| < 1.1 and 
a track with at least one SMT hit, as well as pt > 10 
GeV for Typc-1, and p T > 5 GeV for Typc-2 tau lep- 
tons. There are also requirements concerning overlaps of 
objects: ATZ(^i,t) > 0.4 and A7^(e,r) > 0.4, where r, [i 
and e are as defined above, except that muons that pass 
the overlap criterion do not have to pass the additional 
isolation requirement. To distinguish decays from jets, 
we use a neural network discriminant (T3 |. NNh, and to 
distinguish Type-2 t/j from electrons, we use an addi- 
tional neural network, NN e . We require NNh > 0.9 for 
T\ and T2, and NN C > 0.2 for T2. 

To model the multijet contribution to final states with 
Th decays, we select events with candidates as above, 
but with 0.3 < NN h < 0.8. 



Muons 



E. Jets 



Muons are identified in the muon system, and then 
matched to tracks. They are required to have |?7| < 1.5 
and pt > 15 GeV. The track requirements include a se- 
lection on DCA < 0.02 (0.2) cm for tracks with (without) 
hits in the SMT, where DCA is the distance of closest ap- 
proach of the track to the PV in the transverse plane. 

We require muons to be isolated, meaning that the 
sum of the transverse energies in calorimeter cells in an 
annular region (0.1 < AH < 0.4) around the muon track, 
and the sum of the tracks pr in a cone of A7Z < 0.5 
around the muon track must both be less than 2.5 GeV. 

To estimate the multijet background in the single muon 
sample, we use control samples where the isolation vari- 
ables are required to be between 2.5 GeV and 8 GeV. All 
other criteria are the same as in the signal data sample. 

Because the muon p? is estimated by the pr of the 
matching track in the central tracker, the momentum 
resolution decreases with increasing pt- To restrict the 
analysis to muons with well measured momenta, we 
require the significance of its pr measurement to be 
(1/pt)/o'(1/pt) > 3, where <t(1/pt) is the uncertainty 
on the measurement of the track curvature (inverse of 
the muon track's pt)- This effectively limits muons to 
p T < 200 GeV. 



D. Tau Leptons 

Tau leptons can decay to ev e v Tl fiiy^v T , or hadrons hv T 
(t/j). It is difficult to determine whether a light lepton 
in an event originated from a r, but the signature from 
Th — > hv T differs significantly from that of a jet. The 



We reconstruct jets within \r]\ < 2.5, using an iterative 
midpoint cone algorithm [14| with cone radius of 0.5 and 
a minimum pt requirement of 20 GeV after applying a jet 
energy scale ( JES) correction as discussed in Sec. IIV B 31 
Jets separated from a or an electron by A7Z < 0.5 are 
removed from consideration. 



F. b-jets 

Bottom and charm quarks can travel measurable dis- 
tances from the PV before decaying, so that their decay 
products originate from an identifiable secondary vertex. 
This provides a way of tagging jets coming from a 6(c)- 
quark decay by examining the associated tracks [l5[ . Be- 
fore applying any 6-tagging criteria, the jets are required 
to pass both calorimeter criteria outlined in Sec. IIII El 
and the taggability criteria. A jet is taggable if it is 
matched to a track jet, which is a jet formed from tracks, 
reconstructed using a simple cone-clustering algorithm of 
A1Z < 0.5. At least two tracks are required, with at least 
one having pt > 1 GeV and another with pt > 0.5 GeV. 
Every track in the jet is required to have at least one hit 
in the SMT detector, a DCA < 0.2 cm, and a distance 
of closest approach along the z axis of < 0.4 cm. 

All taggable jets are subjected to a neural network b 
tagging algorithm (l5j whose input variables include the 
DCA of each track in a jet and information on secondary 
vertices in the jet. We define fo-jet candidates by requiring 
that the neural network output be greater than 0.775. 
This algorithm selects about 60% of b jets with pt = 50 
GeV, and only 1% of light flavor (u, d, s quarks or gluon) 
jets. 
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G. Missing Transverse Energy 

Neutrinos or other weakly-interacting neutral particles 
do not leave energy deposits in the detector. Their pres- 
ence is inferred from the measurement of significant $t in 
the event. The missing transverse energy is determined 
from energies deposited in all calorimeter cells. The $t 
is corrected for JES, measured muon px, electron and Th 
energy scales. The JES corrected $t vector is obtained 
by adding the difference between the vector sums of un- 
corrected and JES corrected jet momenta to the uncor- 
rected JpT vector. The muon correction reflects the fact 
that muons deposit little energy in the calorimeter, and 
adjusts the $t for the px of the muon. Finally, electron 
and Th energy corrections arc applied to the appropriate 
calorimeter cells in the J^t calculation. 



IV. MODELING SM PREDICTIONS 
A. SM Event Generation 

We generally estimate SM processes with MC gener- 
ated events. A model-independent search incorporates 
many different processes to properly model the data. We 
use two generators for this purpose, ALPGEN [lj[ for gen- 
eration of all processes, except for diboson production 
which is generated with pythia [13] ■ pythia is also used 
for hadronization and showering. 

ALPGEN uses exact matrix elements at leading orders 
for QCD and electroweak interactions. The benefit of 
using ALPGEN comes from its ability to calculate exact 
leading order terms for processes that include high jet 
multiplicities. ALPGEN produces parton-level events with 
information on color and flavor, and can be matched to 
PYTHIA for parton evolution and hadronization. 

Matching of a parton from ALPGEN to PYTHIA show- 
ering has the fundamental difficulty of separation of the 
hard interaction from initial-state radiation (ISR) and 
final-state radiation (FSR). To address this problem we 
use the MLM matching scheme [lH . In this scheme each 
final state parton from the matrix element is matched in 
A7?. to an evolved jet. We further reject events which 
contain an additional jet not matched to a final state 
parton, except in the sample with the highest number of 
final state partons. 

The following processes are considered, where j is a 
light jet (g,u,d, or s), £ is a lepton, N is an integer > 
and Ip represents a light parton: 

1. W + Nj 

2. Zh* + Nj 

3. W + cc + Nj 

4. W + bb + Nj 

5. Z/ 7 * +cc + Nj 



6. Z/j* + bb + Nj 

7. ii^ {2l + 2v + 2b) + Nj 

8. it -> {lv + 26 + 2lp) + Nj 

9. WW 

10. WZ 

11. ZZ 

The processes involving heavy flavor (HF) quarks (c 
and b) are treated separately from light quark processes 
because they are often associated with particularly in- 
teresting final states, and we generate large number of 
MC events for these final states. Some of these processes 
are included in the light parton simulations, so we remove 
the events with heavy flavor quarks from the light-parton 
samples so as to avoid double-counting. 

For some objects, other programs provide more accu- 
rate simulations of their properties and decays. Specifi- 
cally, TAUOLA [lj| is used for r decays, and evtgen [20[ 
is used for the decay of b hadrons. 

We assume a mass of 172.5 GeV for the top quark, 
consistent with recent measurements [2lj . 

B. Detector Simulation 

The events produced from the above combination of 
generators are processed through the DO detector simu- 
lation and combined with random beam crossing events 
taken from data (Sec. IIVB lj) . The detector simulation 
is based on GEANT 3.2.1 [22j, to which two types of cor- 
rection factors are applied. The first type of correction is 
event reweighting, where an overall correction is applied 
to the MC event, rather than to the measured kinematic 
properties of reconstructed objects. For example, we ap- 
ply weights to account for the difference in reconstruction 
efficiencies between data and MC. Another type of cor- 
rection modifies the objects in a MC event to account for 
the fact that the simulation has better resolution and a 
different energy scale than the detector. These correc- 
tions generally depend on properties of the objects in an 
event. The specific corrections used in this analysis are 
described below. 



1. Instantaneous Luminosity Reweighting 

The trigger selecting random beam crossings records 
data with a different instantaneous luminosity profile 
from that of the triggers utilized to record the data used 
in this search. A weight is introduced in the MC events to 
match the instantaneous luminosity distribution in data. 



7 



2. Zpv Reweightmg 

Our simulated events have a narrower zpy distribution 
than is observed in data. We therefore apply a weight to 
each event, based on the zpy of the event, to increase the 
relative weight of events farther from the center of our 
detector to match the observed distribution. 



3. JES 

We apply JES corrections to jets in both data and MC 
[23| . The purpose of the JES corrections is to correct the 
measured jet energy to that of the particles in the jet. Jet 
energies initially determined from the calorimeter cell en- 
ergies do not exactly correspond to the energies of hnal 
state particles that traverse the calorimeter. As a result, 
a detailed calibration is applied separately in data and 
MC. In general, the energy of all final state particles in- 
side the jet cone, Ej tcl , can be related to the energy mea- 
sured inside the jet cone, E h by Ef A = (E 3 - 0)/(R S). 
Here, O denotes an offset energy, primarily from addi- 
tional interactions in or out of time with an event. R is 
the average response of the calorimeter to the particles 
in a jet, and S is the correction factor for the net energy 
loss from particles that scatter out of or into the jet cone. 
For a given cone radius, O and S are functions of the jet 
rj within the detector. O is also a function of the number 
of reconstructed event vertices and the instantaneous lu- 
minosity; R is the largest correction factor and reflects 
the lower response of the calorimeter to charged hadrons 
relative to electrons and photons. It also includes the 
effect of particle energy loss in front of the calorimeter. 
The primary response correction is derived from studies 
of 7+jct events, and depends on jet energy and pseudo- 
rapidity. For all jets that contain non-isolated muons, we 
add the muon momenta to that of the jet. Under the as- 
sumption that these muons are from semileptonic decays 
of b quarks, we also add an estimated average neutrino 
momentum assumed to be collinear with the jet direction. 

4- Jet Shifting, Smearing, and Removal (JSSR) 

Additional corrections beyond the JES are needed to 
take into account threshold and resolution effects for jets. 
The JSSR corrections are determined from Z/7 — > ee 
+ 1 jet events. The Z/j and the jet should be pro- 
duced approximately back-to-back in (f> with the same 
Pt- This is quantified by a pt imbalance variable, 
AS = (pip — p^^^j /Pt 1 ■ For jets with a px well above 
the reconstruction threshold, the distribution of AS is 
Gaussian in both data and MC. The difference in the 
means of these distributions yields a shift that is applied 
to the MC jet energies to match the data, and a smear- 
ing is applied to MC jets based on the difference in the 
standard deviations of these distributions. Jets that fail 



the pt > 20 GeV requirement after shifting and smearing 
corrections are removed from further consideration. 



5. Efficiencies 

The efficiency of the MC simulation of our detector 
tends to be larger than the true efficiency of the detector. 
To account for this, we introduce scale factors to adjust 
the MC efficiency to match that observed in data. The 
efficiencies for electrons and muons are obtained using 
Z — > ee and Z — > events. One of the decay products 
of the Z boson is the tag object, which is required to pass 
restrictive reconstruction requirements and be matched 
to an object that could have fired the trigger for the 
event. Object efficiencies are then obtained using the 
second object from the Z decay. 



6. Track Pt Resolution 

Electron energies are measured in the calorimeter. 
However, energy deposition does not depend on the 
charge of the electron, which is determined by the cur- 
vature of the associated track in the magnetic field. An 
incorrectly reconstructed track can therefore lead to an 
incorrect charge assignment. Brcmsstrahlung from elec- 
trons can affect the curvature of the tracks. Also, a soft 
interaction in the inner detector can result in the process 
e + — > e + e~e + , leading to charge misidentification if the 
wrong sign electron track is associated with the electron. 
This difficulty is also present in tau decays when at least 
one hadron is produced. 

Because the rate of charge misidentification is not 
properly modeled in the detector simulation, we add a 
scale factor to electron and tau MC events to approx- 
imate the appropriate rate of charge mis-identification. 
We determine this scale factor by using dielectron events 
consistent with Z — > ee decays; and we only consider 
events with dielectron invariant mass between 70 to 110 
GeV to avoid biases against physics beyond the SM. 

The disagreement in track resolution between the data 
and MC also affects muonpT measurement, which is cor- 
rected using smearing parameters determined by com- 
paring the data and MC mass peaks for Z — > fifi and 
J /ip — > decays. 



7. Electron Energy Smearing 

In the simulation, the electron pt reconstructed in the 
calorimeter has a better resolution than in the data. We 
correct this using a Gaussian smearing function tuned to 
reproduce the shape of the Z — > ee peak. 
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8. Jet Taggability 

The jet taggability rates (Sec. IIII F| are found to be 
different for MC and data. To correct for this difference, 
correction factors arc applied as scale factors depending 
on pt, ?7 and zpy of the jet [24j |, 



9. b-taggmg Rate 

As detailed in Sec. IIII F\ we apply a tagging algorithm 
to both data and MC jets to select jets originating from 
heavy (b/c) quarks. However, the algorithm can select 
mistagged light jets. The tagging rates (for both heavy 
and light parton jets) depend on the pt and r\ of the jets. 
The heavy-quark tagging rates are measured separately 
in both data and MC using dedicated samples. The per- 
formance of the 6-tagging algorithm in MC events is bet- 
ter than in data. To correct the tagging rates in MC 
events, we first determine the flavor of the tagged jet by 
matching it in AT? with the initial parton. Depending on 
the flavor of the jet, we apply a per-jet scale factor given 
by SF = e data (p T ,v)/t MC (PT,v), where e data (p T ,?7) and 
e MC (pT,v) are tne ^-tagging efficiencies for a given par- 
ton flavor for data (MC) events. To maintain correct 
normalization, a small downward correction is applied to 
non-&-tagged jets. 



10. Weak Gauge Boson pr 



seven states (e + jets, /i + jets, ee, pp, fie, er, /it ) were 
each selected to correspond to a specific SM process. 

TABLE I: Inclusive final states and their ob- 
ject selections, where p™ m is the minimum al- 
lowed value of pt and |7y| max is the maximum 
allowed value of \rj\. 



Final State Object pf in (GeV) 


1 „ 1 max 

\v\ 


- jets + X a 


e 


35 


1 1 
1.1 




jet 


20 


2.5 






20 




- jets + X b 


[1 


25 


1.5 




jet 


20 


2.5 




E T 


20 




ee + X c 


e 


20 


1.1 


m + x d 


/'• 


15 


1.5 


\xe + X e 


/'• 


15 


1.5 




e 


15 


1.1 


er + X* 


e 


15 


1.1 




T 


15 


1.1 


\1T + X s 


/'• 


15 


1.5 




T 


15 


1.1 



X / e, fi, t, 7 

b X ^ e, fj,, t, 7 

C X y^U, T, 7 

d X jt e, r, 7 
e X^r, 7 
fX^y 
9 X ^ e, 7 



The p T distribution of the Z boson from ALPGEN MC 
is corrected to match the distribution observed in data 
in Z — > ee decays [25[ . A modified reweighting is carried 
over to the W boson px based on the theoretical ratio of 
the W to Z pt spectra [H| . 

11. Acj> 

We apply a A0-dcpcndent weight specifically for this 
analysis using the inclusive distributions described in 
Sec. [V] to correct the A0 between leptons in dilepton 
final states and the lepton and $t in singlc-lepton + jets 
final states. This reweighting affects not only the A<fi dis- 
tributions, but also other quantities that depend on the 
angular distribution of particles such as the pt of the W 
boson. 



V. INCLUSIVE FINAL STATES 

The seven inclusive non-overlapping final states are 
specified in Table U by the relevant objects and their se- 
lection criteria. The additional objects (X in the table) 
are selected as shown in Table |TTJ Events with photons 
are rejected, mainly due to difficulties in modeling. The 



TABLE II: Criteria required for inclusion as 
additional objects (X) in one of the seven fi- 
nal states listed in Table [T] 



Object p 1 ^ 


(GeV) 


^max 


e 


15 


1.1 


l< 


15 


1.5 


T 


15 


1.1 


jet 


20 


2.5 



• e + jets 

The electron + jets final states have more back- 
ground from multijet events, where a jet is misiden- 
tificd as an electron, than the other electron final 
states. Therefore the likelihood criterion used is 
tighter than in other final states, C e > 0.95. We 
also require at least one jet having Et > 20 GeV, 
E T > 20 GeV, and an e p T > 35 GeV. This fi- 
nal state is dominated by W + jets events with 
W — > ev decays. 

• fi + jets 

The pi + jets final state is dominated by W + 
jets events with W — >• pv decays. To reduce the 
amount of multijet background, at least one jet 
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having Ex > 20 GeV is required, as well as Ex> 20 
GeV and a muon with px > 25 GeV. Just as the e 
+ jets final state, this final state is inclusive in jets 
with no other additional objects allowed. 

• ee 

The dielectron final state requires each electron to 
have pt > 20 GeV and C e > 0.85. The electrons 
are also restricted to be in the central calorimeter, 
|?7| < 1.1, and the jets have the same criteria as for 
the other final states. This final state is dominated 
by Z/j* — > ee events. 

• W 

The dimuon final state requires at least two muons 
with the muon-px criteria lowered to pr > 15 GeV 
because of the smaller contribution from multijet 
background. Any jet must have px > 20 GeV. This 
final state is inclusive in both jets and muons, but 
an additional e or r lepton places the event in the 
fie or \it final states. Analogous to the ee channel, 
this final state is dominated by Z/j* — > fifi events. 

• fie 

The fie final state is inclusive except for r leptons; 
e/xr events are assigned to the er final state. This 
final state is dominated by Z/j* — > tt events. 

• er 

The er sample is inclusive in all objects. The elec- 
tron and Th Pt are required to be at least 15 GeV. 
The electron likelihood is set to L e > 0.95 to reduce 
the large multijet background as many apparent Th 
correspond to misidentified jets. The parameter 
that separates electron from hadronic taus, NN e , is 
set to 0.8 to reduce the contribution from dielec- 
tron events. This final state is also dominated by 
Z/"f* — ¥ tt events. 

• fiT 

The fiT state contains at least one muon and one 
Th- It is inclusive in all objects except electrons, 
whose presence would move the event to the er 
final state. This final state is also dominated by 
Z/j* — ¥ tt events. 

VI. INCLUSIVE NORMALIZATION FITS 

Our model does not provide proper normalization 
of different MC contributions because, for example, of 
higher-order corrections needed for the leading-order or 
leading-logarithm cross section calculations. To avoid 
uncertainties in normalization, we perform a fit, de- 
scribed below, for each of the inclusive final states to 
obtain scale factors that reproduce the distributions of 
the selected data using a combination of the SM MC 
and multijet predictions determined from data. We treat 



the Drell-Yan (D-Y) contributions to the ee and fifi fi- 
nal states without light partons separately from those 
with light partons because it improves agreement be- 
tween data and MC. 

The fits for normalization factors arc performed on 
kinematic distributions of different object quantities, al- 
tering the overall normalization of each input process 
contributing to final state so that the x 2 probability 
for that final state is minimized for the combined fit. 
To avoid fitting to data at the highest values of px, 
where new physical processes can be important, we only 
use events that are not in the high px tail, which con- 
tains 10% of the events. Distributions of basic quan- 
tities such as Ex, Pt, V, A</>(obj, Et) of leptons and 
jets (here obj refers to the momentum of the object con- 
sidered) are used in the fits while more complex vari- 
ables are used to check the quality of the overall fit. 
The latter variables include the mass or transverse mass 
AIt = \J{pt,i +Pt,2) 2 — (J>t,i + Pt,2) 2 of two or more 
objects, jet multiplicities, and the px of the W and Z 
bosons. If an event contains any object outside the px 
range defined above, then none of the objects in the event 
are used in the fit. 

The list of the seven final states, the processes that are 
normalized through the inclusive fits to each of the final 
states, and the number of events in each final state are 
shown in Table IIII1 Once the fitted values are extracted, 
the distributions are rescaled accordingly, and the total 
background contribution, B, for a particular final state 
is 

Nbkg 

i 

where the scale factor (Si) for each background process 
(Bi) is determined from the final state in which its con- 
tribution is most important and that scale factor is used 
in all other final states to which that background con- 
tributes. Nbkg refers the total number of all the SM pro- 
cesses contributing to a particular final state. 

A simplified example for the e + jets + X final state 
(X 7^ e, fi, t, 7) is used to illustrate the procedure. The 
e + jots + X state is dominated by W — > ev events, 
but there is a significant contribution from multijet and 
Drell-Yan events. We use the normalization factor for 
the Drell-Yan process, determined through a separate fit 
to the ee + X final state (X 7^ fi, r, 7), in the e + jets 
fit. We also fix the scale factors to one for rare processes 
which have contributions that are too small to fit accu- 
rately in e + jets, such as the tt contribution. We then 
fit for the SM W boson and multijet contributions in the 
data. The fit optimizes agreement between the distri- 
butions in data and the SM prediction for the variables 
listed above. The result of the fit is two overall weights, 
one for W — > ev and one for multijet — > e + jets. 

The distributions of the variables for the input pro- 
cesses are not varied, only their relative contributions. 
The fit is performed using the MINUIT program [27j . For 
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singlc-lcpton states and hadronic r final states, multijct 
events are a significant background. We assume that the 
contribution from other SM processes modeled by the 
MC samples to the multijet background is small. The 
scale factors of input processes for the MC events should 
also account for the contributions of the processes to the 
multijct background. The main effects of contributions 
from any of the MC processes to the multijet background 
would be to decrease the scale factor for backgrounds 
modeled by MC. 

The main purpose of the normalization process is 
to assure that the fundamental SM processes are well- 
modeled. The results of the fit are then checked for qual- 
itative agreement with the data. The overall scale factors 
are checked to compare to those from dedicated analyses. 
If the normalization factors are properly included in the 
MC, then all the scale factors should equal unity. One 
histogram that is included in the overall fit and one check 
histogram that is not part of the fit are shown for each 
of the seven final states in Figs. |TJ [T] In the figures, 
the leading and second electron are the electrons with 
highest pt in the event and next highest pt in the event, 
with a similar definition for leading and second muons 
and jets. 

The electron px distribution in Fig. [TJ shows a clear 
disagreement between data and simulation in this kine- 
matic region arising from the need for a large multijet 
contribution at low pt, and other variables that provide 
better agreement with a smaller multijct contribution. 
However, the discrepancy at low px should not mask the 
presence of new physics at high px, which is the main 
focus of this analysis. 



VII. EXCLUSIVE FINAL STATES 

After determining the normalization scale factors, the 
seven inclusive subsets are merged to create an input file 
for the vista algorithm Each MC and background 
event is given a weight calculated from the data based 
scale factors and any required corrections. The vista al- 
gorithm, developed by the CDF Collaboration, is a tool 
that performs a broad check of the agreement between 
data and the SM. We modified the CDF algorithm for 
our analysis strategy as described above. The resultant 
VISTA@D0 algorithm focuses on the DO high pt data to 
determine whether the data can be adequately described 
by the SM or if significant discrepancies can be con- 
firmed, vista mainly examines discrepancies that affect 
the overall distributions rather than narrow regions of 
phase space, addressing the numbers of expected events 
and MC / data agreement across full distributions of cho- 
sen variables. 

The use of standard object identification criteria 
(Sec. IIII[) provides great simplification in the analysis as 
data can be partitioned into exclusive final states. The 
events are separated into homogeneous subsets of events 
according to the objects contained in each event, rcsult- 
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FIG. 1: (color online) e + jets final state (a) electron pr 
histogram and (b) transverse mass {e,$r) check histogram. 



ing in 117 exclusive final states. Examples of such exclu- 
sive final states include /i ± r T + 2 jets + $t , e fi* + 2 
jets + $t j e + e+ + 3 jets, and /x + 4 jets + $t- 

VISTA performs two types of checks: first, it does a 
normalization-only check on the number of events in 
each exclusive state; the goodness of the fit is calcu- 
lated using Poisson probabilities. Second, it calculates a 
Kolmogorov-Smirnov statistic (and resulting fit probabil- 
ity) for the consistency of all the kinematic distributions 
in any final state with the predicted SM distributions. 
Both of these results require additional interpretation be- 
cause of the large number of trials (number of final states 
and/or the number of distributions) involved. When ob- 
serving many final states, some disagreement is expected 
from statistical fluctuations in the data. Thus the Pois- 
son probability used to determine agreement is corrected 
to reflect this multiple testing. A similar effect occurs 
when comparing kinematic distributions, and again the 
probabilities are first converted to standard deviations 
and then corrected for the number of distributions ex- 
amined. 

Another algorithm we use to search for new physics is 
called sleuth Q , developed at the DO experiment during 
Run I (1992-1996) of the Tevatron. sleuth is an attempt 
to systematically search for new physics as an excess at 
the largest values of J^Pt- This variable corresponds 
to the sum of the values of the scalar px of all objects 
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FIG. 2: (color online) p + jets final state (a) $t histogram 
and (b) transverse mass (p, $t) check histogram. 



in the event, including the Tfir- The SLEUTH algorithm 
is quasi-model independent, where "quasi" refers to the 
assumption that the physics beyond SM will appear as 
an excess of events at large pr- Therefore sleuth is ex- 
pected to be most sensitive to high-mass objects decaying 
into relatively few final-state particles. 

For SLEUTH, the VISTA exclusive and 1-jet final states 
are merged, as are the 2 and 3-jet final states, and light- 
lcpton universality is assumed, combining eX and [iX 
channels. Making these assumptions greatly reduces the 
number of states considered in sleuth relative to vista, 
and thus the trials factor, improving the statistical sensi- 
tivity by diminishing the chance of observing a large fluc- 
tuation. Next, the ^pt distribution in each channel is 
scanned to find a cutoff that maximizes the significance of 
any excess in data relative to the SM background, defin- 
ing a lower bound for the ^2pt selection. Finally, the 
probability for consistency with the SM of the largest 
values of the px is corrected for the number of possi- 
ble lower bounds in any distribution, and subsequently 
for the number of final states examined by sleuth. This 
corrected probability corresponds to the probability that 
any individual final state would yield probabilities as 
small as observed. We define a significant output from 
SLEUTH as one with a corrected probability of < 0.001 
(that is over 3 Gaussian standard deviations from the SM 
prediction using a one-sided confidence interval). 



FIG. 3: (color online) ee final state (a) leading electron (with 
highest pr) Pt fit histogram and (b) invariant mass (e,e) check 
histogram. 



VIII. SENSITIVITY TEST 



To check the sensitivity of a search with sleuth, we 
examine whether a top quark (produced in it pairs) which 
contributes objects with high pt would have been discov- 
ered in the current data sample. For this test, we used all 
the background samples, except for the it MC. The main 
concern is whether other final states would compensate 
for the missing it events, and thus SLEUTH would not be 
sensitive to it production in data. 

We examine the £jjbb$T final state, which we expect 
to be dominated by it events. Figure [5] shows that pres- 
ence or absence of a it signal has a great impact. With a 
threshold of 0.001, the sleuth test, including the it MC, 
yields a statistical probability of compatibility of 0.98 af- 
ter correcting for the number of trials. However, without 
the it contribution this probability is < 1.1 x 10~ 5 . In 
Fig. [8]and other sleuth plots, the insets show the results 
for data and MC that pass the ^ px cut maximizing the 
significance of excess in data. 
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FIG. 4: (color online) fiji final state (a) leading muon pr fit 
histogram and (b) invariant mass (p,, p) check histogram. 



FIG. 5: (color online) pe final state (a) electron pr fit his- 
togram and (b) invariant mass (fi, e) check histogram. 



IX. RESULTS 

A. Numerical discrepancy using the VISTA analysis 

In vista, the separation of the input data into final 
states completely defined by the objects in an event, 
yields a total of 117 unique exclusive final states. The 
probability (V) that the yield observed in data results 
from a statistical fluctuation of the SM sample in chan- 
nel f s is determined from 



V = l-(l-p fl 



f sa N fs x Pfa 



(2) 



where Nf s is the number of trials and pf s is the probabil- 
ity that the number of events predicted for the channel 
f a in the SM would fluctuate to what is observed in data, 
before applying the correction for the number of trials. 
The number of trials is Nf a =117, corresponding to the 
number of final states, and 



Pfa 
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(N - N B f 



2a% 



dN 



OO 

E 

Wdata 



-N 



(3) 



where Ng and <tb are the expected SM event yield from 
background and its uncertainty, respectively, and iVdata 
is the number of events observed in any channel. The 



Gaussian significance is the value of <j that satisfies the 
equation 



1 



2tt 



" dx 



(4) 



The final state probabilities converted into standard de- 
viations, before the correction factor for the number of 
trials, are shown in Fig. |H1 This distribution shows most 
final states near a = 0, with some excess for a > 3. Of 
the 117 final states, two show significant discrepancy af- 
ter correction for the number of trials. These are the 
final states p, + 2 jets + $t, with a probability corre- 
sponding to a 4.5 a discrepancy, and ii + [i~ + $t with a 
discrepancy of 6.7tr (also shown in Fig [9]). 

The discrepancy for the n + 2 jets + JpT final state 
shows the greatest difference from the SM prediction in 
the modeling of jet distributions. There is a significant 
excess in the number of jets at high \rj\, which points to 
likely problems with modeling ISR/FSR jets in the for- 
ward region, as can be seen in F ig. [10k . This difference is 
observed in dedicated analyses [281 ] , and the discrepancy 
becomes less severe when using SHERPA [29[ MC events. 

The fi + p~ + $t discrepancy can be attributed to dif- 
ficulties modeling the muon momentum distribution for 
high pt muons. As noted in Scc. lIVB6l the muon smear- 
ing modeling is based on muons from Z and J/^p decays, 
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FIG. 6: (color online) er final state (a) The Ac/>(e,^T) fit 
histogram and (b) invariant mass (e,r) check histogram. 



FIG. 7: (color online) /xt final state (a) muon pr fit histogram 
and (b) invariant mass (fi, r) check histogram. 



dominated by muons below 60 GeV, and is not as reliable 
at high pr . The prime signature of poorly simulated high 
Pt muons is an excess of $r because of the mismodcling 
of the resolution of the mismeasured track. The A</> be- 
tween the positive muon and $t in the + $t final 
state is shown in Fig. [TUb . where the excess tends to be 
for events where the $t is collinear with a muon. 



B. VISTA Shape Analysis of Discrepancies in 
Distributions 



3<7 level after correcting for the trials. The majority of 
these are related to spatial distributions involving jets. 
All these discrepancies are related to known simplifica- 
tions in our modeling assumptions, e.g., no systematic 
uncertainties taken into account, aside from the adjust- 
ments made by the normalization factors. These dis- 
crepancies would not be expected to severely affect the 
sleuth search for new physics at high pr tails. All 16 
discrepant distributions are shown in Figs. H2H15l and arc 
listed in Table ITVl In the figures, the second jet refers to 
the lower pr jet in the two jet final states. 



The 117 final states contribute a total of 5543 indi- 
vidual one-dimensional distributions in various variables, 
and comparison between simulation and data is per- 
formed for each. The trials-factor adjusted probability 
is determined from V = 1 — (1 — p s h P ) 5543 , where p s h p 
is the Kolmogorov-Smirnov (KS) probability to observe 
a discrepancy for any individual distribution (before ap- 
plying the correction for 5543 trials). As with the prob- 
ability for a final state normalization discrepancy in any 
final state, the probability for a discrepancy in a spec- 
trum is converted into units of standard deviation. Any 
deviation >3cr is considered discrepant. The distribution 
of deviations before correction for the number of trials is 
shown in Fig. [TT1 

Sixteen distributions are found to be discrepant at the 



(J. SLEUTH 

All vista final states are used as input to SLEUTH, and 
the 117 inclusive final states are folded into 31 final states 
after applying global charge conjugation invariance, re- 
binning in the number of jets, and assuming light lepton 
universality. The two vista final states that show broad 
numerical excesses are found again with the sleuth al- 
gorithm, as expected. No additional final states have a 
significant sleuth output, as defined in Sec. IVIII 

In the sleuth runs performed at CDF, using a slightly 
different analysis strategy, the four most interesting ob- 
served final states were fJ^e , ^ ± e ± + 2 jets + ^r, /J, e 
+ pr, and £±1^ l' + $ T in 2.0 fb" 1 of integrated lumi- 
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TABLE III: The contributions used in the in- 
clusive fits for each of the inclusive final states 
and the number of data events in each final 
state. The dominant SM process is listed first 
for each final state. In the er and fir final 
states, the multijet background also includes 
a contribution from W + jets. 



State 



SM process 



Events 



e + jets + X 



jj, + jets + X 



ee 



fie 



(IT 



X 



X 



X 



X 



X 



W + jets 
Multijet 

W/Z + HF 

W + jets 
Multijet 

W/Z + HF 

D-Y + Olp 
D-Y l-31p ratio 

W/Z + HF 



D-Y 
D-Y 



lp 
1-3 lp 



W/Z + HF 

Z — y tt 
Multijet 

W/Z + HF 

Z — y tt 
Multijet 

W/Z + HF 

Z — y TT 

Multijet 
W/Z + HF 



40k 



50k 



25k 



24k 



0.34k 



1.3k 



1.0k 



nosity. These states were also among the most discrepant 
observed by CDF in 0.9 pb _1 [|| of integrated luminosity. 
Our results for these states are shown in Figs. [T51[T71 and 
IT51 except for /i ± e ± + 2 jets + $t, for which we find no 
events with 0.16 events expected. Figure [T9l shows the 
similar final state, where the muon and electron are of 
opposite sign rather than of the same sign where CDF 
sees a discrepancy. None of these states are significantly 
discrepant in our analysis. 

The sleuth final states with V < 0.99 are shown in 
Table [V] A plot including all of the final state probabil- 
ities converted to units of a can be seen in Fig. The 
final state + r T + $t, which was not identified as hav- 
ing a significant discrepancy between data and the SM 
expectation in vista, falls close to our sleuth thresh- 
old. Figure |2"T1 shows the J^Pt distribution for this final 
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FIG. 8: (color online) Sensitivity to new physics test using the 
tt final state, (a) The tt MC is included, yielding only minor 
differences between data and SM background. The statisti- 
cal agreement between the data and MC for the distribution 
shown on inset is nearly 2 a. (b) The results of the entire 
analysis without the it MC. In this case, sleuth passes the 
criterion of interest at 0.001 for this final state. The insets 
shows the distribution beyond the ^2pt cutoff. "Other" refers 
to contributions too small to list, including W + bb — > eubb 
events, W + cc^t luce events, W + lp — > tv + lp events, and 
diboson events. 



state. 



X. CONCLUSIONS 

We have done a global study of DO high pr data to 
search for significant deviations from the standard model 
expectations. This broad search for BSM physics has 
been done on 1.1 fb _1 of integrated luminosity collected 
in Run II of the Fermilab Tevatron Collider in the DO ex- 
periment. Using the vista algorithm, a total of 117 ex- 
clusive final states and 5543 kinematic distributions were 
compared to the SM background predictions. Only two 
out of 117 exclusive final states, fj, + 2 jets + $t an d 
+^t, show a statistically significant discrepancy. 
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FIG. 9: (color online) Distribution of discrepancies for the 
117 final states relative to the SM in terms of standard de- 
viations calculated in VISTA final state before accounting for 
the trials factors. The curve represents a Gaussian distribu- 
tion centered at zero to guide the eye. The distribution is 
expected to obey Poisson statistics, which is the reason the 
distribution is narrower than the Gaussian. 



FIG. 11: (color online) The a distribution for the 5543 VISTA 
comparisons before accounting for number if the trials. The 
curve represents a Gaussian distribution centered at zero to 
guide the eye. There are 116 distributions in the underflow 
bin with a < —10. This is expected as histograms with KS 
probabilities > 0.99999 are rounded to 1, and appear in the 
underflow bin. 
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FIG. 10: (color online) (a) The distribution of the pseudora- 
pidity of the second jet with respect ot the center of the DO 
detector in the /i + 2 jets + $t channel. "Other" contains 
distributions too small to list individually, W + bb, diboson, 
tt, and D-Y + Ip. (b) The A(j> distribution between the 
and the $t for the +-^t final state. "Other" contains 

distributions too small to list individually, diboson and tt. 



TABLE IV: The full list of VISTA results with dis- 
crepant distributions listed by final state. 

vista Final State Histogram a 



2 jets + IS T 



M T {W,:h) 4.4 

ATZ(fx,j 2 ) 4.4 

M(jij 2 ) 4.0 

Arj(ii,i 2 ) 3.8 



1 jet + $ T 



p T (W) 8.1 
SpT 5.1 
Pt(u) 4.1 
Mt^^t) 4.1 
A<Mm,j) 3.1 



e ± + 2 jets 


+# T Ar?(ii,i 2 ) 4.2 




M T (j 2 ,$ T ) 4.0 




M T {W,j 2 ) 3.0 


e ± + 1 jet - 


-#r A<t>(e+,j) 5.5 




p T ( e ±) 4.4 




p T (W) 3.8 




IfjT 3.1 



Given the known modeling difficulties in both final states 
together with our neglect in this study of systematical 
uncertainties, we cannot attribute the observed discrep- 
ancies to sources of physics beyond the standard model. 
A quasi-model-independent search for new physics was 
also performed using the algorithm sleuth by looking 
for statistically significant excess at high "YjPt in a wide 
array of exclusive final states. No additional final states 
cross the discovery threshold in sleuth beyond the ex- 
cesses noted by vista. 
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TABLE V: The sleuth states with V < 0.99. The 
value of V represents the corresponding probability 
without taking into account the trial factor. 



Final State 


7> 


-pa 


£+£- + $ T 


< icr & 


< 0.001 


+2j + # T 


< 10~ 5 


< 0.001 




3.9 x 1CT 5 


0.0050 


e ± + p T + ij 


0.00036 


0.019 


eV T +26 + $T 


0.0028 


0.12 


^T± + 2j + ^ T 


0.0028 


0.12 


£± + 26 + #r 


0.0077 


0.3 


e*^ + ^ T 


0.0081 


0.31 


£ ± T ± 


0.057 


0.91 


£ ± + 26 + 2j + # T 


0.099 


0.98 



' The value of V is not necessarily accurate below 0.001. 
The important check is whether the value drops below 
the threshold. Further discussion can be found in [g|. 
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FIG. 12: (color online) The discrepant distributions in the n 
+ 2 jets + ]pT exclusive final state, (a) The transverse mass 
distribution of the W boson plus second jet, (b) the A1Z be- 
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FIG. 13: (color online) The discrepant distributions in the fj, + 1 jet + $t exclusive final state: (a) the pr of the W boson, 
(b) the sum of the scalar values of pr of the fi, jet, and |?t, (c) the transverse mass of the fj, and $t, (d) the pt of the /x, and 
(e) the A(f> between the muon and the jet. "Other" contains distributions too small to list individually, W + bb, diboson, tt, 
and D-Y + Ip. 
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FIG. 14: (color online) The discrepant histograms in the e + 
2 jets + $t exclusive final state, (a) The A77 between the 
two jets, (b) the transverse mass of the trailing jet and $t, 
and (c) the transverse mass distribution of the W boson plus 
trailing jet. "Other" contains distributions too small to list 
individually, diboson, D-Y, and ti 
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FIG. 15: (color online) The discrepant distributions in the e 
+ 1 jet + $t exclusive final state, (a) The A<j> between the e 
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FIG. 16: (color online) Check of most discrepant CDF plots 
from H, /i ± e i 
Epr cut. 



The inset shows the distribution above the 



FIG. 19: (color online) Since there are no data events in the 
fj, e + 2 jets + $t final state, the distribution for fj, + 
2 jets + TfiT is shown. The inset shows the distribution above 
the EpT cut. "Other" contains the Z — > /j/j and W/Z + bb 
distributions. 
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FIG. 17: (color online) Check of most discrepant CDF plots 
from [s|, fj, ± e ± + $t- The inset shows the distribution above 
the EpT cut. 



FIG. 20: (color online) Distribution of final state sleuth 
probabilities converted into units of a before inclusion of the 
final state trials factor. 
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FIG. 18: (color online) Check of most discrepant CDF plots 
from + $t- The inset shows the distribution above 

the EpT cut. 
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FIG. 21: (color online) sleuth plot for p- + t t + $ T . The 
inset shows the distribution above the Epr cut. "Other" in- 
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